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I. INTRODUCTION 

Duriig the last ten years the Wilson-Kadanoff ap¬ 
proach to the Renormalization Group (RG), based 
on the block spin concept, has been the subject of a re¬ 
vival in both Statistical Physics and Field Theory. This 
originates in recent developments iii which have 
now turned it into an efficient tool, the so-called effective 
average action method 0], allowing to investigate non¬ 
perturbative phenomena. This method implements on 
the effective action T - the Gibbs free energy - the idea 
of integration of high-energy modes that underlies any 
RG approach. The whole method consists in building 
an effective free energy Tfe at scale k for the high-energy 
modes that have been integrated out and in following its 
evolution with the scale k through an exact equation Q . 
The main drawback of this equation is that it cannot 
be handled in actual calculations without truncations of 
Tfe. It is thus of utmost importance to know whether 
the truncations used provide converged and accurate re¬ 
sults. As it is well known, the problem of convergence is 
also crucial in perturbation theory where it requires re¬ 
summation of series. Let us emphasize that this problem 
is far from being solved in general since Borel summabil- 
ity, which is the key point to resum perturbative series, is 
not TOnerically proven and may even turn out not to hold 
(see la [3 for a review). It is then important to dispose of 
an alternative method, not relying on an expansion in a 
coupling constant and thus not requiring a priori resum¬ 
mation. Good indications on the convergence properties 
of the effective average action method have been already 


provided by its ability to tackle with highly nontrivial 
problems, such as low-energy QGD 0, the abelian Higgs 
model relevant for superconductivity the p hase di¬ 
agram of Hes M, frustrated magnets |1 lL hd m, 
the Gross-Neveu model in three dimensions [id , the 
randomly dilute Ising model IIaI , the Kosterlitz-Thouless 
transition , etc (see |2n| for a review and for 

an exhaustive bibliography). A systematic investigation 
of the convergence and accuracy issues is however still 
lacking. 

We propose here, on the example of the three- 
dimensional Ising model, to study the convergence and 
optimization of the accuracy of the effective average ac¬ 
tion method truncated both in derivatives, at order d^, 
and in fields. We study, in particular, the role of the cut¬ 
off function used to separate the low- and high- energy 
modes, on the determination of the critical exponents v 
and 77 . 

In section II, we briefly introduce the basic ideas un¬ 
derlying the effective average action method. We then 
discuss in section III the truncations necessary to deal 
with concrete calculations. We motivate, in section IV, 
the use of the Principle of Minimal Sensitivity (PMS) 
to optimize the results. Then, we apply this technique 
successively within the Local Potential Approximation 
(LPA), section V, and at order of the derivative ex¬ 
pansion, section VI. 


II. THE EFFECTIVE AVERAGE ACTION 
METHOD 
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Historically, the block spin concept was first imple¬ 
mented, in the continuum, on the Hamiltonian. This 
procedure consists in separating, within the partition 
function, the microscopic fields into a high- and a low- 
energy part and in integrating out the high-energy part 
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to get an effective Hamiltonian for the remaining low- 
energy modes. The iteration of this procedure gener¬ 
ates a sequence, a flow, of scale-dependent Hamiltonians, 
parametrized by a running scale k, and describing the 
same long distance physics. The critical properties are 
then determined by the behavior of the system around 
the fixed point of the flow of Hamiltonians @]. However, 
due to technical difficulties this nonpertur- 

bative renormalization procedure has been mainly used 
as a conceptual basis for perturbative calculations rather 
than as a practical tool to investigate nonperturbative 
aspects of field theories and critical phenomena. The sit¬ 
uation has ch ange d when it has been reali zed, ma inly by 
Ellwanger [2l|2a|M|pL MorrisSMIllllili^ 
and Wetterich |3,HISllli that, rather than 
the Hamiltonian H, one should consider the effective ac¬ 
tion r - the Gibbs free energy - as the central quantity 
to be renormalized. In the spirit of the original Wilso¬ 
nian formulation one builds a running effective action 
Tfe that only includes high-energy fluctuations with mo¬ 
menta > k'^. This implies that, on the one hand, at 
the underlying microscopic scale fc = A, T^ coincides with 
the classical Hamiltonian H since no fluctuation has yet 
been taken into account. On the other hand, when the 
running scale is lowered to fc = 0, i.e. when all fluctu¬ 
ations have been integrated out, the standard effective 
action T is recovered. To summarize, T^ continuously in¬ 
terpolates between the microscopic Hamiltonian H and 
the free energy: 


r Tfc^A = H 

( 1 ) 

{ = r . 

Since, by definition, Tj, is built up from the high-energy 
fluctuations of the microscopic system, the low-energy 
modes - with < k^ - must be removed from the 
running partition function. This is most easily achieved 
by adding to the original Hamiltonian a scale-dependent 
mass term A.Hk. Then, the running partition function 
with a source term writes [20): 

Zk[J] = J e~^^x] - ^Hklx] + J-X (2) 
with J.x = f d'^q J{q)x{-q) and 

^Hk[x] = \ j Rk{q)x{q)x{-q) (3) 

where x(g) is the microscopic field. In Eq.®, Rk{q) is 
chosen in such a way that it acts as a cutoff function 
that decouples the low- and high- energy modes. This 
imposes several constraints: 


Rk{q) ' 

- k^ for q^ < k^ 

(4) 

Rkiq) - 

0 for . 

(5) 


Equation 0 means that, at low-momentum with respect 
to k, Rk{q) essentially acts as a mass, i.e. an infrared 
cutoff, which prevents the propagation of the low-energy 
modes. This ensures that these modes do not contribute 
to Tfc [^. Eq.® implies that Rk{q) does not affect the 
propagation of high-energy modes. They are thus almost 
fully taken into account in Zk and, consequently, in T^. 

In order to recover the limits (1), Rtiq) must also sat¬ 
isfy: 


Rk(q) —> oo when k ^ A at fixed q (6) 

which ensures that Tfe coincides with the microscopic 
Hamiltonian H when fc —> A, and 

Rk{q) 0 identically, when fc ^ 0 (7) 


which ensures that, in the limit of vanishing k, one re¬ 
covers the standard effective action T. Note that since 
we are only interested here in the universal long distance 
behavior and not in quantities depending on microscopic 
details, we send A to oo. 

The effective average action T/; is then defined as: 


Tfc [(/)] = - In Zfc [ J] + J.(^ - AHk [0] (8) 


where (j) stands for the running order parameter (j)k{q)'. 


<l^k{q) = {xiq))k 


6\nZk[J] 

SJ{q) 


J=Q 


(9) 


It follows from the definition 0 that r/c[( 5 i] essentially 
corresponds to the Legendre transform of In Zk [ J], up to 
the mass term AHk which allows to recover the limits 

m E3- 

The effective average action T^ follows an exact equa¬ 
tion which controls its evolution with the running scale 

fci: 


dtTkm = 11^ d,Rk{q) {rf [</>(<z)] 

( 10 ) 

(‘2'\ 

where t = ln(fc/A) and T). '[(j)] is the second functional 
derivative of Tfc with respect to the field 4‘{q). We em¬ 
phasize that Eq. cn is exact and thus contains all per¬ 
turbative and nonperturbative features of the underlying 
theory (see 113 for technical details and | 23 | for a review 
of the applications of this equation to concrete physical 
issues). 


III. TRUNCATIONS OF THE EFFECTIVE 
AVERAGE ACTION 

Equation cni) is a functional partial integro-differential 
equation that has obviously no known solution in the 
general case. Therefore, to render it tractable, one has 
to truncate the effective action T^. The most natural 
truncation, well suited to the study of the long distance 
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physics, is the derivative expansion. It consists in writing 
an ansatz for as a power series in dcj)- Let us first 
consider the case of an 0{N) invariant theory for which 
the ansatz at the order writes : 

TfcM = J i Uk (p) + ]^Zk{p) + 

^ . ( 11 ) 

Wkip) {d^pf + 0 {d^)\ 

where </> is a TV-component vector and p = is the 

0{N) invariant. In Eo. llTTl) . Uk (p) corresponds to the 
potential part of Fj, while Zk (p) and Yk (p) correspond to 
the field renormalization functions. Thus, with Zfe(p) = 1 
and Ffc(p) = 0 , Ea. dllll provides the ansatz for the so- 
called Local Potential Approximation (LPA) where the 
anomalous dimension vanishes. This kind of ansatz has 
been successfully used in several cases among which the 
0{N) and Gross-Neveu models HE HI However, 
to deal with more complicated models, e.g. with matrix¬ 
like order parameters, a further approximation is almost 
unavoidable [llL IlI 13^ . Indeed, when the symmetry is 
lower than 0{N), there are several invariants and the 
number of independent functions analogous to Zk{p-, ■■■) 
and Yk{p, ...) grows. In this case, the integration of the 
flow can be very demanding. It is then very convenient 
to further truncate the functions C/fc (p,...), Zk{p,...) in 
power series of p and of all other invariants. 

Here, we focus on the Ising model, described by a 
scalar, Z2-invariant field theory, considered as a toy 
model to study the derivative and field expansions. In 
this case, since the only independent field renormaliza¬ 
tion function is Zk{p), the function Yk{p) can be set to 
zero. The field truncation then writes: 

n 

Uk{p) =^C/,,fc(p-poy 

( 12 ) 

P 

Zk{p) =y^^Zi^k{p - poT 

V 2 = 0 

where po = 4 >o being a particular configuration of 

the field (p. We shall come back on this point later. 

The truncation in fields conveys two nice properties. 
Eirst, with the ansatz CH and the RG flow equa¬ 
tion m leads to a finite set of ordinary coupled differen¬ 
tial equations for the coupling constants Ui^s and Zi^s 
simpler to solve than the partial differential equations 
obeyed by the full functions Uk{p) and Zk{p). Second, 
even the lowest order approximations, in which only the 
first nontrivial terms of Uk (p) and Zk (p) are kept, give a 
fairly good qualitative picture of the physics |33| • 
However, the study of the truncated version of Eq. ns 
raises several important questions: 

i) Does the derivative expansion converge and does it 
provide a satisfying accuracy at low orders ? The ques¬ 
tion of the convergence of the derivative expansion, in 


its full generality, has not yet been considered and ap¬ 
pears to be a major and open challenge. In practice, 
one is less interested in this delicate question than in 
the quality of the results and their improvement as the 
order of the derivative expansion is increased. In the 
case of 0 {N) models, very accurate results have been ob¬ 
tained at second order in the derivative expansion. Eor 
instance, Wetterich et al. have shown that handling the 
full field-dependence of the potential Uk{p) and of the 
field renormalization functions Zk{p) and Yk[p) leads to 
results that can compete with the world best estimates, 
at least for the critical exponent v [^. The value ob¬ 
tained for the anomalous dimension p is less accurate. Its 
definition being linked to the momentum dependence of 
the two-point correlation function, an accurate determi¬ 
nation of T] probably requires higher order terms in the 
derivative expansion. This question will be investigated 
in a forthcoming article |^. 

ii) Does the field expansion of Uk{p) and Zk{p) con¬ 
verge and how rapidly ? Once again, the general question 
of convergence has not yet been investigated. Neverthe¬ 
less, several works have dealt with field truncations at 
high order within the LPA ^ 3 ; HE 1 ^ or with 
a field independent field-renormalization | 3 ^. i.e. with 
Zk[p) = Zo^k- They suggest that a few orders suffice 
to obtain reasonably converged values of critical expo¬ 
nents. To our knowledge, their computation using also 
an expansion of Zk (p) has been only studied in the Ising 
model and using a power-law cutoff function . In this 
study, we extend this analysis to two other families of 
cutoff functions, leading to more accurate results. 

The questions i) and ii) are linked with a corollary 
issue, which resides in the choice of cutoff function. One 
naturally inquires about its influence, and in particular: 

Hi) Gan the accuracy be improved through the choice 
of cutoff function Rk ? Of course, when no truncation 
is made, an exact solution for F[0] = limfe^o Lfe [</>], does 
not depend on the function Rk used, whereas any kind 
of truncation induces a spurious dependence on it. One 
can thus wonder how to optimize the choice of this cutoff 
function. This question is not as trivial as it seems since 
one has to decide of an optimization criterion: rapidity 
of convergence of the expansions in powers of derivatives, 
fields or amplitudes HE HE HE ? accuracy of the 
results ? sensitivity of the results with respect to the 
cutoff ? We specifically concentrate on these two latter 
issues in the following. 


IV. OPTIMIZATION AND PRINCIPLE OF 
MINIMAL SENSITIVITY 

Up to now, attempts to optimize nonperturbative RG 
equations have been mainly worked out in the Polchin- 
ski equation [d^], in particular at second order in the 
derivative exp ansion. For instance, Ball et al. and 
Gomellas have tried to suppress the cutoff and nor¬ 
malization dependence of the exponents ly and rj by us- 
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ing the Principal of Minimal Sensitivity (PMS) [^. We 
shall not pursue within this framework since it has now 
been widely recognized that the effective average action 
method is the most efficient way to deal with the nonper- 
turbative RG. We will thus consider this latter formalism. 

In the context of the effective average action method, 
within the framework of LPA, Litim h as prop osed to con¬ 
sider the quantity C, defined by QH EE 1^ EE ■ 


rnin ( rf^[(/-(g)] 

9 >0 


4>—<Pq 


+ Rkiq) = Ck 


(13) 


(21 

where P ).' [(j){q)]+Rk{q) is the inverse of the full regular¬ 
ized propagator and C parametrizes the gap amplitude. 
According to Litim, the gap is bounded from above and 
the best cutoff functions are those which maximize this 

gap lillilliilililEE: 


C'opt = niax(C') when varying Rk- (14) 

The idea behind this criterion is that the larger C, the 
more stable the truncated RG flow. Indeed, it has been 
shown that the maximum of the gap corresponds to the 
largest radius of convergence of an amplitude expansion. 
This suggests that the optimal selected regulators should 
have nice properties such as improvin g th e convergence 
of the field expansion [43. l43l lift 1^ . Moreover in 

[iE it has been shown that, within the LPA, the criterion 
da is also linked to a PMS. 

At this stage, let us shed the light on some impor¬ 
tant features of the “gap criterion”. First Eq. Cl does 
not select a unique cutoff function: many Rk maximizing 
the gap have been exhibited for instance in Also, 
the various optimized cutoff functions, solutions of Eq. 
d, can lead to quantitatively different critical expo¬ 
nents depending on the specific properties of a given cut¬ 
off function, like its asymptotic behavior (see below and 
compare 0 and 01 )■ The quality of the results there¬ 
fore relies on the choice of the type of optimized regula¬ 
tor. Second, beyond the LPA, the implementation of the 
gap criterion appears to be nontrivial. Indeed, the 
field renormalization function Zk{po) induces an implicit 
i?fc-dependence in P^ l4‘{q)] complicates the maxi¬ 
mization of the gap. Moreover, it is not completely clear 
whether, beyond the LPA, this criterion would still con¬ 
vey the nice properties it shows at the lowest order of the 
derivative expansion and, in particular, its link to a PMS. 
As we are specifically concerned here with the question 
of the sensitivity of the results with respect to the cut¬ 
off function, we favor a method that directly probes the 
dependence of the critical exponents on the cutoff func¬ 
tion. We have decided to base our analysis on the PMS, 
that can always be simply implemented and has already 
proven its efficiency. 

Let us recall how it works. Suppose, for instance, that 
we compute a quantity Q in an approximate way. The 
approximation used may induce a dependence of Q on 
a parameter - denoted here a - which is spurious. The 


PMS consists in choosing for a the value a pms for which 
Q is stationary: 


dQ{a) 


da 


= 0 


(15) 


One thus expects that imposing such a constraint, sat¬ 
isfied by Q computed without approximation, improves 
the approximate determination of this quantity. The ob¬ 
vious drawback of this method is that Eq. 03) can have 
many solutions. This worsens if several quantities are si¬ 
multaneously studied, and lead to distinct solutions. An 
additional criterion is then necessary to select a unique 
one. 

We first study the LPA of the scalar, 22-invariant field 
theory relevant for the Ising model. We show that the 
PMS allows one to optimize the quality of the results. 
We then study the O(d^) approximation of the deriva¬ 
tive expansion and show that the PMS leads to accurate 
results provided we add some new inputs to discriminate 
the solutions. 


V. THE LOCAL POTENTIAL 
APPROXIMATION OF THE ISING MODEL 


Let us recall that the LPA consists in approximating 
Tfe by: 


Tk[(t>] = j |[/fc(p) -f i (16) 

i.e. in neglecting the field renormalization. This ansatz, 
once plugged into Eq. enables to get the evolution 
equation for Uk- Actually, working with dimensionless 
quantities is necessary to get a fixed point, so that we 
define: 


/ ^ Rkiq^) q^ 

Uk = k-'^Uk 
, p = k'^~‘^p . 


The RG equation obeyed by Uk writes: 


duk 

~dt 


duk + {d- 2)pu'k - VdLQ{w) 


(17) 


(18) 


where Uk = Uk{p), = 2‘^+^7r'^/^r((i/2), prime means 
derivation with respect to p^ w = u'f. + 2pu'l, and 





‘2y'^r'{y) 

y[l + r{y)] + w' 


(19) 


The nonperturbative features of the evolution of the po¬ 
tential are entirely encoded in the integral Lg, called 
threshold function [^ . 
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We now study Eq. within a field truncation: 

n 

Uk{p) =^u.,{p-pQy ( 20 ) 

i=l 

where we have suppressed the index k for the coupling 
constants. Once Uk{p) is truncated at a finite order n 
of the field expansion, the field configuration po around 
which it is expanded matters. Two configurations have 
been widely studied: the vanishing field configuration, 
Po = 0, and the configuration where Uk{p) has a nontriv¬ 
ial minimum |^: 


duk 

dp 



( 21 ) 


All the studies performed using field truncations show 
that the convergence properties are improved by expand¬ 
ing around the minimum rather than around the zero 
field configuration Therefore, we choose the for¬ 

mer. 

We also need to choose families of cutoff functions Rk 
to perform calculations. For simplicity, we restrict for 
now our study to families of cutoff functions Rk depend¬ 
ing on a single parameter. We extend this to a two- 
parameter family in section IVI PI We consider two usual 
cutoff functions. The first one is the exponential cutoff, 
which has been often used and constitutes an efficient 
and robust regulator |^. The other one, the theta cut¬ 
off, has been introduced by Litim ^5. It presents the 
advantage of leading to threshold functions that can be 
analytically computed. We extend theses functions, by 
multiplyiM them by a factor a, to two one-parameter 
families [^Id^: 


fexp,a{y) ~ ® gy _ ^ 
r0,aiy) ^ ~ 


( 22 ) 


Note that both of these cutoff functions can be opti¬ 
mized according to the gap criterion. 

For both families, we investigate the a dependence of 
the critical exponent v over a large range of a, for each 
order n of the field expansion, up to the tenth power of 
p. We indeed expect the most relevant operators to be 
contained in the first terms, and thus the evolution of v as 
a function of the order of the truncation to be stabilized 
at, or before, the tenth order. 

We find that, at each order, v{a) exhibits a single point 
of minimal sensitivity for both cutoff functions. For rg a, 
(FIG.HJ, the minimum occurs at apMS = as in [d^ . 
with an optimized v equal to v{apMs) = vpMS = 0.650. 
For Texp.a, (FIG.IH, one has apMS = 6.03 and vpMS = 
0.651 (see Table I). Both cutoff functions lead to very 
similar optimal results for p, differing by less than 0.5% 
to all orders n, as shown in FIG.|21 The converged values 



0.2 0.6 1 1.4 1.8 

a 



FIG. 1: Curves i^(a) for the cutoff rg^c,, for different trunca¬ 
tions of the potential Uk{p)- Note that, for n > 6 - lower 
figure -, the v axis is magnified. 


of V are reached below the percent level in both cases after 
only a few orders (n = d), as expected. 

The insert of FIG.|2| where the evolution of p with the 
order n of the field expansion is compared for a = 0.1 
- chosen for illustrative purpose - and a = apMS = 
shows that the same convergence level is reached inde¬ 
pendently of a right from the n = d order, though the 
asymptotic values of j^(O.l) and v{apMs) differ signifi¬ 
cantly. This shows that the rapidity of convergence cri¬ 
terion is helpless here to select a cutoff. 

We now compare our results with those obtained 
through the gap criterion. As displayed in FIG. Q 
OiPMS = 1 exactly with rg^a to all orders. For this cut¬ 
off function the op m.s value coincides with that given 
by the gap criterion |^. For rexp.a, apMS converges to 
6.03 (see lower curves in FIG. 0 ) whereas the gap crite¬ 
rion selects an optimal parameter oq = 3.92 |^. In this 
case, the two methods seem to differ. However, since 
the variations of v, once converged, do not exceed one 
percent in the whole range a S [oi ~ 1.2,02 — 7d], 
we do not expect the two methods of optimization to 
lead to drastically different critical exponents. Indeed, 
p(ao) = v{oiPMs) up to 10“"^. Thus for the two families 
of regulators considered here, the PMS and gap criteria 
coincide. It has been argued that this property holds, 
within the LPA, for more general families of regulators 

0143. 
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a 

0.652 
0.6518 
0.6516 
_ 0.6514 
0.6512 
^ 0.651 

0.6508 
0.6506 
0.6504 

23456789 10 

a 

FIG. 2: Curves u{a) for the cutoff rexp.a, for different trun¬ 
cations of the potential Uk{p). Note that, for n > 6 - lower 
figure the v axis is magnified. 

Note that for the exponential cutoff, the standard 
choice a = 1 leads to v = 0.658. This value, which does 
not correspond to an optimized one, differs by a little bit 
more than one percent from vpMS- For completeness, we 
also mention that the power-law regulator optimized via 
the gap criterion — r{y) = y~‘^ — leads to a less accurate 
result: = 0.660 0,^3 • 

Finally, let us emphasize that the world best value v = 
0.6304(13) (see Table D) lies below all curves v{a) for 
both cutoff functions and that the PMS solutions for p 
are minima. Thus, ^{apMs) is the most accurate value 
achievable within each family of cutoff functions studied 
here. The PMS therefore constitutes a powerful method 
to optimize the cutoff function in order to reach the best 
accuracy on the critical exponents. 

VI. ORDER (9^ OF THE DERIVATIVE 
EXPANSION 

We now show how the PMS can be consistently im¬ 
plemented at the order of the derivative expansion 
for which, as far as we know, no optimization procedure 
has ever been implemented within the effective average 
action method. We dispose of two physical quantities 
candidates for a PMS analysis: v and f]. We perform 
both analyses independently, with each cutoff function. 



FIG. 3: vpMS for re^a and as a function of the order 

n of the field truncation of Uk{p). The two curves almost 
superimpose for all n. In the insert, p is displayed for rg^a, 
for both apMS and a = 0.1. 



Ref. 

V 

V 


a) 

0.651 

0 

LPA 

b),b’) 

0.650 

0 


c) 

0.660 

0 


d) 

0.6307 

0.0467 


a) 

0.6281 

0.0443 


b) 

0.6260 

0.0470 


c) 

0.6175 

0.0542 

7-loop 

e) 

0.6304(13) 

0.0335(25) 

MC 

f) 

0.6297(5) 

0.0362(8) 

Exp. 

g) 

0.636(31) 

0.045(11) 


h) 

0.6298(90) 



TABLE I: Critical exponents of the three dimensional Ising 
model, a), b), b’), c) and d) are computed from the effective 
average action method: a) with rexp.apMs (present work); b) 
with rg^apMs (present work); b’) with re,c«=i c) with 

power-law cutoff li^. l53l : d) with rexp,a=i without field ex¬ 
pansion 1^ : e) from perturbation theory including 7-loop 
contributions Ih^ : f) from Monte Carlo simulations IH^ : g) 
from experiment in mixing transition 1^ : h) from experiment 
in liquid-vapor transition (computed from 3p — 2 — a IH^ l. 


We show in section Em that the PMS allows one to 
improve the accuracy on both exponents. We especially 
highlight that accuracy is not synonymous of rapidity of 
convergence of the field expansion. In section rvrm we 
bring out a necessary condition for the independent im¬ 
plementation of the two PMS on p and rj to be consistent. 
We then check that our results meet this condition. In 
section we exhibit cases where, contrary to what 

occurs in the LPA, multiple PMS solutions exist. We 
show that a unique one can be selected thanks to general 
arguments. We end up by extending the analysis to a 
two-parameter family of cutoff functions. 
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A. Accuracy of the PMS solution and convergence 
of the field expansion 

This section is devoted to showing that the PMS is 
still, at order d^, the appropriate tool to find, within a 
class of cutoff functions, the one giving the best accu¬ 
racy. Though it seems counter-intuitive, we emphasize 
that this cutoff function does not coincide with the one 
providing the fastest convergence of the field expansion 
of Zk{p). To this purpose, we implement both PMS in¬ 
dependently on v and 77 , postponing the coherence of this 
to the next section. 

Working with a nontrivial held renormalization func¬ 
tion Zk {p) , dimensionless and renormalized quantities are 
necessary in order to get a hxed point, so that we dehne: 


(23) 


The RG equation 


= pzk + pz'i,{d - 2 + rj) + Vd{z',. + 2pz'^)Lf{w, Zk,v)- 
4vdpzk{3uk + 2pu'l')L^{w,Zk,r])- 
2vd (2 -f 1 /d) p {z'kf Zk, 1]) + 

(4/d) VdP {iu'k + 2pUkYMf{w, Zk,'n) + 

( 8 /d) Vdpz'k{‘iul + 2pUk)M'l'^'^(w,Zk,ri) + 

(4/d) Vd p {z'kfM^+‘^{w, Zk,v) 

(24) 

where w = u'+2j5u", prime means derivative with respect 
to p, and the threshold functions are dehned, for n > 1 , 
by: 


r{y) = 


Rk{q^ 


Zo.k 

UkiP) = k~'^ Ukip) 
p = Zo,fe P~‘^p. 

Zk{p) 


with 


Zk{p) = 




where Z^^k is dehned in Eq. 03- 
obeyed by Zk (p) writes : 


rdr ^ f d d/2-i^y r'{y) +vyr{y) 

Ln(w,Zk,v)=n dyy/ 

do {P{y) + w}^+^ 

do [Piy) + wp 


7/(1 -k r{y) + yp{y)){pr{y) + 2yr'{y)) 


(25) 


P{y) + w 


-2r]{r{y) + yr'{y)) - 4y{2r'{y) + yr"{y)) 
where 

P{y) = y{zk+r{y)). 


(26) 

(27) 


The anomalous dimension rj is given by: 

d=-^lnZo.fe- (28) 

As previously, we truncate the held renormalization 
function Zk {p) up to the p-th power of p: 

p 

Zk[p) =^z,(p-po)b (29) 

i=0 

We use, for the potential Uk{p), the expansion given in 
Ea. lEHIl . up to the term, which represents a very ac¬ 
curate approximation of Uk (p) in the vicinity of its mini¬ 
mum as shown in the previous section. We expand Zk{p) 
up to the ninth power of p which turns out to be sufficient 
to obtain converged results. 
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FIG. 4: Curves i/(a) for Texp.o, for different truncations of 
the field renormalization Zk(p}- For p > 5 - lower hgure -, 
the V axis is magnihed. Note that the curve uio^s shows two 
extrema [62]. 

At each order p of the held expansion of Zk (p), we have 
computed the exponents v and rj as functions of a for 
both cutoff functions and rexp.a. FIG. 0 and FIG. 
inigather the curves representing these functions, labelled 
uwZp, p — 0,..., 9, on the example of rexp.a- They are 
displayed on a range of a around the extremum and sepa¬ 
rated in two distinct hgures since the p > 5 curves would 
be superimposed without magnihcation. This seems to 
indicate that the held expansion converges, at least on 


























a 



a 

FIG. 5: Curves »?(«) for rexp,a, for different truncations of the 
field renormalization Zk(p). For p > 5 - lower figure the p 
axis is magnified. 

the whole range of a studied. The same conclusion holds 
for rg^a, with very similar curves (that we therefore do 
not show), up to the important subtlety, discussed in 
section iwn that two PMS solutions exist in this case 
and that only one has to be considered. We call 
and 7]°^ (a) the two limit functions obtained for p —> oo. 
In practice, we approximate these functions by those at 
p = 9. 

Let us emphasize that for both cutoff functions, i) the 
rapidity of convergence to and 77 °“ (a) and ii) the 

asymptotic values and 77 °“, depend both on a. One 
can thus naturally wonder whether the values of a for 
which the convergence is the fastest coincide with those 
for which the exponents are the most accurate compared 
with the world best results. We shall show that this is 
not the case contrary to what is widely believed. 

t. Accuracy 

We first bring out that the PMS exponents are, as in 
the LPA case, the most accurate ones. We have deter¬ 
mined, for each p, the values cxpMs(p) '^pms(p) 
which, respectively, v and 77 reach their extremum. The 
corresponding exponents are referred to, in the follow¬ 
ing, as = via'^pMsip)) and rfpj^^g = vi^PMsiP))- 

The obtained PMS asymptotic values are — 0-6281 


and rf^MS ~ 0.0443 for the exponential cutoff, and 
^PMS = 0.6260 and ~ 0.0470 for the theta cutoff. 

These values of the exponents are indeed the best achiev¬ 
able within each class of cutoff functions studied, since 
the world best value of v lies above the sets of curves in 
FIG. 21 and since the extremum is a maximum (and vice 
versa for 77 ) (see Tableland HJ). The PMS is thus, 
as in the LPA case, the appropriate tool to find, among 
a family of cutoff functions, the one providing the best 
accuracy. 


2. Rapidity of convergence 

The evolution of v^p^g and Ppj^g with the order p of 
the field expansion of Zk{p) is displayed in FIG.Elfor both 
cutoff functions. The convergence of vpMS and ppMSj at 
the percent level, requires at least p = 4 for both cutoff 
functions. However, there exist values of the parameter 
a, for instance a = 1.80 for rg^a, for which the conver¬ 
gence is faster than for apMS- This is illustrated in the 
inserts of FIG. El Indeed p{a = 1.8) has already con¬ 
verged at the percent level for p = 3, but to a different 
value than Thus, the PMS exponents, which are 

the most accurate, are not those converging the fastest. 

We conclude that i) the PMS leads to the most ac¬ 
curate exponents within each class of cutoffs studied, ii) 
a criterion based on rapidity of convergence of the field 
expansion would be here misleading since it would select 
cutoff functions leading to exponents significantly differ¬ 
ing from the PMS ones. 


B. Consistency condition for independent PMS 
implementations 

We have implemented and discussed the PMS analy¬ 
ses independently on v and p along the previous section. 
This has naturally led us to two distinct PMS values of 
a at each order p, cx'^.^g (p) and cc),]yjg (p). One can thus 
wonder whether it makes sense to compute two different 
quantities with two different cutoff functions. We now 
provide a natural condition for the whole procedure to 
be consistent. 

Let us notice that since the field expansion seems to 
converge (as shown in the previous section), the two 
sequences a'pMsiP) '^pms(p) converge. The 
asymptotic value a^Ms(P ~ (^esp. cxpMsiP ~ °°)) 
is the one that achieves the minimum dependence of the 
exponent 77 (resp. u) on the cutoff function at order of 
the derivative expansion. There is no reason for them to 
coincide. However, the discrepancy between the apMs’s 
does not matter as long as choosing one or the other does 
not change significantly the value of each exponent. A 
consistency condition is thus: 

^(«PMs(oo)) J7(a^Ms(oo)) (30) 
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P 


FIG. 6 : upMS and rjpMS for both rexp.c and re,a as functions 
of the order p of the field expansion of Zk{p)- In the inserts 
are displayed, for re,a, and for two distinct values of a, a = 
otpMs(p) and a = 1.8, the critical exponents v (upper insert) 
and rj (lower insert) as functions of p. 


and 


viapMsi^)) - »?(apMs(oo)) ■ (31) 

Reciprocally, large discrepancies between the values, at 
the two apMs’Sj of an exponent would be an indication of 
a failure of convergence. It could be imputed to either a 
too low order of expansion, or to an unappropriate choice 
of cutoff functions family. 

In principle, we should check the consistency over the 
whole set of exponents describing the model. Let us how¬ 
ever show that once this condition is satisfied by two in¬ 
dependent exponents, it is automatically by all the oth¬ 
ers, provided the scaling relations hold within the cho¬ 
sen truncation scheme (in fields and derivatives). Let us 
first emphasize that it has been observed in all instances 
where it has been studied that the scaling relations re¬ 
main precisely verified order by order in the field expan¬ 
sion, although the exponents vary much with the order. 



P 



P 


FIG. 7: Plot of and ol'pms s-® functions of the order p 

of the field expansion of Zk{p) for both cutoff functions Cexp.a 
and re,a- ctpMS (resp- ctpMS ) i® the value where lies the 
PMS extremum of p, (resp. u). 


We thus assume that computing the critical exponents 
either directly or from the scaling relations is (almost) 
equivalent. In this case, an exponent, 7 for instance, re¬ 
lated to p and 7 through the scaling relation: 


7 (a) = p{a){2 - r]{a)) 
obviously verifies for all a: 


d'y 

da 


dv , 

;r(2-^) 

da 


dr] 

'^d^' 


(32) 


(33) 


In the simple case where apj^^g coincides with a^p^^g, 
we deduce from Eq. that 7 ( 0 ) also reaches its ex¬ 
tremum for this apMS- Thus, a^p^^g = apj^^g = apj^^g 
and the consistency is trivially verified for 7 also. In 
the general case where the apMs’s are distinct, if they 
correspond to consistent exponents p and t] according to 
Ea. lTTnil and ED, one is ensured that both exponents are 
almost stationary between these two apMs’s, provided 
the functions pia) and r]{a) are smooth enough in this 
range. Hence, it follows from Eq. ESJ that dj/da almost 
vanishes both at a = a^j^^g and at a = apj^g. This 
means that 7 (a) is also stationary around these points, 
and thus, 7 computed from a PMS analysis should verify: 


7(apMs(oo)) - 7(apMs(oo)) ^ 7(apMs(oo))> (34) 
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i.e. 7 meets the consistency condition. Using the same 
argument for all the other exponents, we deduce that the 
independent implementations of the PMS on all expo¬ 
nents are consistent once they are for two independent 
ones. 

Let us now examine our results. FIG. 0 sketches 
a PM sip) ^'pMsip) functions of the order p of 
the field truncation, for both cutoff functions rexp.a and 
rg^a- Let us set out a few comments. First, the func¬ 
tions O'PM sip) and <apMs(p) converge as expected. On 
the one hand, oJ,Msip) turns out to be very stable, and 
roughly converging as fast as tfpMS- This originates in 
the very peaked shape of the function 7?(a) (lower curves 
of fig.EI). On the other hand, OpMs shows larger os¬ 
cillations, due to the flatness of the function (lower 
curves of FIG. El). It is worth mentioning that since the 
exponents have almost converged at p = 4 (FIG. 01, the 
fluctuations on the corresponding opMS values induce 
negligible variations on them for p > 4. 

Let us now show that the independent analyses of p 
and give consistent results with respect to Eqs. CT 
and m- The asymptotic values are approximated by 
those at p = 9. The consistency condition is trivially ver¬ 
ified for Texp.a since in this case <apMs(oc) — OpMsi'^) 
(see FIG.EI). For rg^a, we find: 

Wio^PMsi^)) - '^(«PMs(oo))l - 

(35) 

\viopMsi°^)) - viopMsio°))\ - 6-10 

which are both negligible. Thus, in this case also, the 
consistency condition is fulfilled. We draw the conclu¬ 
sion that the PMS analyses have selected a unique opti¬ 
mal value for each exponent v and p although the corre¬ 
sponding opMs’s do not coincide. They enable to deduce 
the remaining critical exponents as well. 


C. Discrimination of multiple PMS extrema 

The results discussed in section are associated 
with a particular PMS solution while several ones can ex¬ 
ist, leading to significantly different exponents [^. This 
happens for rg^a, (see FIG. 0. We now expose the gen¬ 
eral arguments we used to discriminate between the dif¬ 
ferent PMS solutions. 

Suppose that the derivative expansion is studied order 
by order without field truncation (or equivalently that 
the field expansion is perfectly converged). If the deriva¬ 
tive expansion converges, the corrections on exponents 
must be smaller and smaller as the order of the expan¬ 
sion is increased, at least at sufficiently large order. On 
the other hand, as the asymptotic value of any observable 
is exact, it must be independent of the cutoff function. 
Thus, for any quantity, all cutoff functions lead to the 
same asymptotic - exact - value, although not at the 
same speed. In practice, the aim is to reach it as fast as 
possible. This means that, at least beyond a certain or¬ 
der, the best cutoff for the derivative expansion is the one 


which leads to the fastest convergence. Note that this is 
not the case for the field expansion where the rapidity of 
convergence does not provide a criterion to discriminate 
between various PMS solutions. 

Of course, this asymptotic value could be reached only 
after large fluctuations occuring at first orders, as in the 
field expansion (see FIG. 0 . However, contrary to this 
case and provided p is not too large, we expect the first 
orders of the derivative expansion to already lead to re¬ 
liable results. Under this hypothesis, we get two natural 
criteria to select a unique PMS solution when several ex¬ 
ist. The first one consists in keeping, for each family 
of cutoff functions, only the PMS solutions that have a 
counterpart in the other(s) family(ies), i.e. that lead to 
(almost) the same critical exponents. This means in our 
case that we keep only the PMS solutions that verify (in 
obvious notations): 


exp Q 

^PMS — ^PMS 

exp g 

VpMS — VPMS 


(36) 


since these exponents are stationary not only inside a 
family of cutoff functions but also from one family to 
the other. The second criterion consists in applying our 
previous hypothesis of rapid convergence already at order 
d^: we assume that no large fluctuation occurs between 
the LPA and approximation. We thus select the PMS 
solution that minimizes, on the exponents, the correction 
of order to the LPA. 

Both criteria allow one to discriminate between the 
two distinct PMS solutions obtained for p and p with 
rg,a (see the curve itiozg in FIG.0. They happen to pair 
for both exponents, at roughly af^Ms — OpMS — and 
~ of MS — According to the second criterion, 
we exclude the second PMS solutions located at opms — 
6.5, which lead for both v and p to much larger deviations 
than the first ones compared with the LPA result: p = 0 
and p = 0.650 (see FIG.0 curve (rg^a) uiq). The first 
criterion leads to the same choice since i) we have checked 
that with rexp,Q; only one PMS solution exists for p (resp. 
for p), and ii) the corresponding exponent is very similar 
to the one at the first PMS solution for p (resp. for p) 
with rg^a, see Table Hland FIG.0 Thus, our two criteria 
to select a unique PMS solution are consistent [^. 


D. Influence of a second parameter 

In the previous sections, we have restricted our anal¬ 
yses to the influence of the parameter a, amplitude of 
the cutoff functions, on the critical exponents. The op¬ 
timized results obtained with the two families of cutoff 
functions are very close together. It is thus natural to 
test the robustness of this result. In this section we in¬ 
vestigate the influence of other deformations of the usual 
cutoff functions focusing on the exponential cutoff. Two 
generalizations of rexp,ct come naturally. They consist 
in changing i) expy exp/3y and ii) expy expyP 
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FIG. 8 : Curves v{a) and r]{a) for both cutoff functions rg^a 
and Texp.a within LPA (labelled uio) and at O(d^) of the 
derivative expansion (labelled 11102 : 9 ), for the maximal trun¬ 
cations of Uk{p) and Zk{p) computed here. The two PMS 
extrema for rg.c are shown for both ly and p. 

mumiii. The deformation i) reveals actually use¬ 
less since it is equivalent to a rescaling of the running 
scale k in Rk which is immaterial. We hence study the 
two-parameter generalization of Texp: 

T'exp,a,f3{y) = OL ^ (37) 

We perform the full PMS analyses of v and p over 
the two-parameter space spanned by a and (3, within 
the LPA and at order of the derivative expansion, for 
the maximal field truncations of Uk{p) and Zk{p) con¬ 
sidered here. We find a unique two-dimensional PMS 
for both exponents, and at both orders. It lies at 
Q^pms — ^PMS ~ 2.25, P'pMs — I^PMS ~ gives 

Vpms = 0.04426 and vpMS = 0.6281 at order 9^. It 
turns out that our prior choice (3 = 1 was very close 
to (3pMS, and thus the a optimization performed in the 
previous section already enabled us to almost reach this 
minimum. The two-parameter PMS exponents thus dif¬ 
fer by less than a tenth of percent from those obtained 
previously (see Table P). 

For illustration purpose, we isolate in FIG. El the be¬ 
havior of the (3 parameter, fixing a to its PMS value de¬ 


termined in section Em It displays the r\{(3) and v{(3) 
functions, for the converged field truncations. Both ex¬ 
ponents exhibit a single PMS solution for (3 very close to 
one {(3pMS = 1.001 in LPA and (3p]^g = Pp^s ~ 0.993 
at order 9^). 






P 


FIG. 9: Curves p(/3) and ^{(3) for rexp,apMS:P’ within LPA 
(labelled uio) and at 0{d^) of the derivative expansion (la¬ 
belled uiozg) for the maximal truncations of Uk(p) and Zk(p) 
considered here. apMS is the value obtained in section Em 

As shown in FIG. El the P dependence of v and rj is 
quite sharp. It raises a natural question: had we fixed 
P far from Ppms to perform the a PMS analysis, what 
would have we obtained ? In other words, would the a 
optimization have suffice to retrieve exponents close to 
the two-parameter PMS ones? To investigate this ques¬ 
tion, we have fixed 9 = 2, which seems from FIG. Elto al¬ 
ter much ? 7 , and determined OpjvfS ^'pms- corre¬ 
sponding exponents are Pp'^g = 0.05573, Vpp^g = 0.6246 
at order 9^. The discrepancy with the two-parameter 
PMS exponents is quite significant for rj, whereas the 
larger exponents - v and the others computed from the 
scaling relations - only undergo a few percent variation. 
This originates in the difference of nature of both expo¬ 
nents. On the one hand, the exponent v is related to the 
behavior of the mass, embodied in the minimum of the ef¬ 
fective potential. The weakness of the sensitivity of v on 
the cutoff function, at order 9^ of the derivative expan¬ 
sion, suggests that the effective potential is already well 
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approximated at this order, and thus provides an accu¬ 
rate determination of close to the exact value. On the 
other hand, rj describes the momentum-dependent part 
of the two-spin correlation function, for which the order 
truncation constitutes a very rough ansatz. Hence, 
the determination of 77 is rather poor at this order and 
improving it probably requires higher derivative orders. 
This is directly reflected in the non-negligible dependence 
of r] on the cutoff function underlined above. 

The conclusion to be drawn from this is that, as pre¬ 
viously, the PMS is the appropriate method to select, 
among a class of cutoff functions, the one that achieves 
the best accuracy, in so far as it minimizes the distance 
to the world best values for both exponents and at both 
orders. Moreover, the PMS reveals itself all the more cru¬ 
cial that the variations with respect to a given parameter 
are large. 

VII. CONCLUSION 

We have implemented the Principle of Minimal Sensi¬ 
tivity to improve critical exponents within the framework 
of the nonperturbative RG. We have shown that it always 
allows to reach the most accurate results achievable in 
the class of cutoff functions under scrutiny. Within the 
LPA, the PMS exponents turn out to almost coincide 
with those obtained through the principle of maximiza¬ 
tion of the gap, and the method is easily generalizable at 
order . 

Two main drawbacks are usually attributed to the im¬ 
plementation of the PMS: i) several solutions of the PMS 
can exist and render its implementation ambiguous, ii) 
it is not clear whether it indeed improves the results. 
We have shown on the example of the Ising model, that 


within the context of the effective average action method, 
these drawbacks either can be circumvented or do not 
exist at all. We have indeed brought out that a unique 
solution of the PMS can always be selected, thanks to 
very reasonable criteria, and furthermore this solution 
represents the most accurate determination of the crit¬ 
ical exponents. The PMS thus appears as a safe and 
powerful method to optimize the results obtained in the 
nonperturbative RG context. An important and rather 
unexpected aspect of our analysis is that the rapidity of 
convergence of the field expansion is not optimal where 
the accuracy is. 

Let us also emphasize that, even within a rather mod¬ 
est truncation involving the potential expansion up to 
order and the field renormalization expansion up to 
order the accuracy reached on v is below the percent 
level compared with the world best results. This suggests 
that, with the same kind of computational complexity, a 
comparable accuracy can be achieved for more compli¬ 
cated models. 

Finally, the determination of 77 is poorer, which is to be 
imputed to the roughness of the ansatz to describe the 
full momentum dependence of the two-spin correlation 
function. Improving it is likely to require inclusion of 
terms of order d‘^. This will be investigated in [^. 
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the relevant eigenvalue of the stability matrix, linearized 
around the fixed point. This computation differs from the 
one implemented in m. It generates a small discrepancy 
(0.5%) between the (d) value in Tabled] and i/(a = 1) = 
0.6272 computed here, see lower curves in FIG.|1] 

[62] Note that several PMS solutions can lead to almost de¬ 
generate critical exponents. See for instance the curve 
uio 25 on FIG. Elfor rexp.o where two PMS exist at Pi and 
P 2 . In this case, whichever point can be selected arbitrar¬ 
ily, since anyway the discrepancy between up ms {Pi) and 
upMs{P2) is negligible (it does not exceed a few tenths 
of percent here). 

[63] Actually, for a —> 00, both u and rj approach an asymp¬ 
totic value for rexp.a that, by extending the notion of 
PMS to infinite a could be considered as a second PMS 
solution. However, the values of both v and rj thus ob¬ 
tained - v{a = c») ~ 0.60 and r]{a = 00) ~ 0.124 - 
are far from those at the second PMS of re,a - iz ~ 0.61 
and r} ~ 0.088 - and therefore cannot be considered as 
consistent. 


